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ABSTRACT 

We have obtained a high spatial resolution X-ray image of the nucleus of NGC 
1068 using the High Resolution Camera (HRC-I) on board the Chandra X-ray 
Observatory, which provides an unprecedented view of the innermost 1 arcsec 
radius region of this galaxy. The HRC image resolves the narrow line region into 
Q_j X-ray emission clumps matching bright emission-line clouds in the HST [OIII] 

A5007 images and allows comparison with sub-arcsec scale radio jet for the first 
time. Two distinct X-ray knots are revealed at 1.3-1.4 arcsec northeast and south- 
west of the nucleus. Based on the combined X-ray, [OIII], and radio continuum 
morphology, we identify the locations of intense radio jet-cloud interaction. The 
[OIII] to soft X-ray ratios show that some of these clouds are strongly affected 
by shock heating, whereas in other locations the jet simply thrusts through with 
no signs of strong interaction. This is further strengthened by the presence of 
a kT ~ 1 keV coUisionally ionized component in the ACIS spectrum of a shock 
heated cloud HST-G. We estimate that the kinematic luminosity of the jet-driven 
shocks is 6 X 10^^ erg s~^, a negligible fraction (10~^) of the estimated total jet 
power. 

Subject headings: X-rays: galaxies — galaxies: Seyfert — galaxies: jets — galax- 
ies: individual (NGC 1068) 



1. Introduction 

NGC 1068 is the brightest prototype Seyfert 2 active gal actic nucleus (AGN), and als o 



one of the closest, with a distance of only 14.4 Mpc (1" = 72 pc: lBland-Hawthorn et al.lll997l ). 
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This makes NGC 1068 an optimal target to investigate the morphology and kinematics of 
the different components in the AGN central region s: the infrared-emitt i ng molecular gas 
the X-ray obscuring torus, an d the ionized gas (e.g. 
Storchi-Bergmann et al.ll2012l ). 



Exposito et al.ll201ll : iKrips et al.ll2011 



Several detailed studies have been performed at multi-wavelength with the highest spa- 
tial resolution instruments avai lable, probing the nuclear region down to scales of the order 
of, or smaller than ~10 pc (e.g.. lEvans et al.lll99ll : ICapetti et al.lll997l : lGallimore et al.ll2004j : 



Groves et al.ll2004f l. 



he picture emergi ng from these studies is that of a complex circum- 



nuclear environment (iRaban et al.ll2009l ). Previous X-ray studies identified a compact pho- 
toionized component, unresolved on scales of ~ 100 pc, emitting the soft X-ray spectrurn 



studied through d ispersive grating spectroscopy ( iKinkhabwala et al.l 120021 : lOgle et al.l 12003 
Evans et al.ll2010( ). and an ionization cone on ~kpc scales, observed in HST [OIII] images 



(Young et al. 


2001; 


Oele et al. 


2003) 



Besides reporting a general go o d cor respondence between the [OIII] distribution and 
the X-ray morphology, lYoung et al.l ( 120011 ) found that the ACIS -S spectra of sev e ral di ffuse 
emission regions could not be fitted with thermal plasma models. iBrinkman et al.l ( 120021 ) and 
Ogle et al.l ( 120031 ) demonstrated that the X-ray emission of the centermost region (unresolved 
with ACIS) and that of a cloud NE of the nucleus is dominated by photoionization. This 
would also explain the general agreement between X-ray and [OIII] morphology. 

While these conclusions are strong, we may be still missing a part of the story. In partic- 
ular, the radial velocities of the emission lines suggest forces at work w ithin &' of the nucleus, 
accel erating outflows in the central 2" and then gradual deceleration (ICrenshaw fc Kraemer 
2000|). These motions suggest clouds accelerated by nuclear winds, then decelerating because 
of interaction with the a mbient ISM ; thes e velocities may a lso reflect the interaction of the 
radio jets with the ISM (lAxon et al.lll998l : lOgle et al.ll2003l ). The inner few hundred parsec 
region is where jet-ISM interactions are expected to be most intensive, but this region has 
not been investigated in great detail in X-rays. Previous deep Chandra ACIS observations, 
while providing illuminating X-ray findings, suffered significant photon pileup in the r < 2" 
(< 144 pc) nuclear region. Even the extended emis sion is affected by pileup in the ACIS 
exposure with a 3.2 s frame time (lYoung et al.ll200ll ). 



In this paper, we report on the X-ray morphology mapped with a resolution of 0.26", 
~ 22 pc, in the nuclear region obtained from the deconvolved HRC imaging of NGC 1068, 
which is free from pileup effect. Taking advantage of archival ACIS image with subpixel 
resolution, we are able to identify X-ray structures with features seen in other wavebands 
and locations of intense jet-cloud interaction. Detailed imaging and spectral analysis of 
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X-ray emission from the nucleus to galaxy scale will be presented in a forthcoming paper, 
utilizing all available Chandra ACIS and HRC data. 



Observations and Data Reduction 



NGC 1068 was observed on October 2, 2010, with the Chandra HRC-I (IMurrav et al. 



20001) for 40.18 ks (ObsID 12705, PI: Fabbiano). The nominal pointing was (a = 02M2'"41."28, 
5 = — 00°00'34."9) with 82.9 degree roll angle. The HRC data were reprocessecQ with the 
Chandra Interactive Analysis of Observations software package (CIAO; version 4.2) tool 
hrc_process_event^ using Chandra Calibration Database (CALDB) version 4.3.1. The 
net exposure time was 39.7 ks after screening for brief periods of elevated background. Since 
there is a known artifact^] in the Chandra/HKC point spread function (PSF), we have used 
the CIAO tool make_psf .asymmetry _region to identify its location as precaution. The ar- 
tifact is in the SE quadrant (position angle [P.A.]~87° to ~137°), not related to any of 
the morphological features in NGC 1068 that we resolve and disc uss in later section s. The 
archival ACIS S3 observation (ObsID 370, previously presented in lYoung et al.ll200ll ) shown 
in this work was taken on February 22, 2002 with 0.1 s frame time and a total exposure 
of 1.43 ks. We reprocessed the ACIS data to generate new le vel 2 file that has the latest 
calibration and the subpixel resolution algorithm ("EDSER". iLi et al.ll2003l ) applied. To 
evaluate the accuracy of astrometry, we compared the positions of X-ray point sources to 
the coordinates frorn the the Naval Observatory Merged Astrometric Dataset (NOMAD; 
Zacharias et al.l 120041 ) and derived absolute astrometric accuracy of 0.3" (Icr). 



3. Imaging and Spectral Analysis 



Following our Chandra HRC study of the Seyfert 1 galaxy NGC 4151 ( jWang et al. 



2009al ) , we performed ima ge restoration using the expectation thr ough Markov Chain Monte 



Carlo (EMC2) algorithm (lEsch et al.ll2004J : iKarovska et al.l 120051 1. The effe ctiveness of this 



method wa s demonstrated with a number of astronomical imaging studies (IKarovska et al. 



2005, 


2007; 


Wang et al. 


2009a) 



The Chandra PSF was simulated with the Chandra Ray 



^ http : //cxc . harvard . edu/ ciao/ th.reads/createL2/ 
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Tracer (ChaRTO) using a monochromatic energy at 1 keV (recommended for HRC data; see 
ChaRT thread noted above). The rays were then projected onto the HRC detector with 
CIAO tool psf _project_ray adopting a 0.2" Gaussian bhirring to generate the PSF image. 

The EMC2 restored (500 iterations) HRC image of the NGC 1068 nuclear region (the 
central 8" x 8") is presented in Figure [H which also shows the raw HRC-I image, and the 
ACIS-S images at native pixel resolution and at subpixel resolution. 



3.1. X-ray Morphology of the Narrow-Line Region 



Although all four Chandra images show similar northeast-southwest elongation, finer 
features are seen in the subpixel rebinned ACTS image, which were not disce rnible in the 
raw A CTS image due to undersampling with the native pixels (c.f. Figure 2 in lYoung et al. 
200ll ). These higher resolution features are also evident in the raw HRC image and more 
strikingly in the PSF-deconvolved HRC image. 

The mos t probable location of the nucleus i s geri erally accepted to be the VLBI radio 
source "SI" (iGallimore et al.l Il996al : iKishimotd Il999l ). We registered our X-ray image so 
that the X-ray peak emission is coincident with UV peak in (jCapetti et al.l 119971 ). and 0.1" 
north of SI (the active nucleus), consistent with an in direct view of the nucleus in the X-rays 
(obscured by a column density exceeding 10 



cm ^, 



Matt et al.ll2000f ). 



The X-ray emission in the northeast (NE) direction, where the optical extended narrow 
line region is, shows a cylindrical morphology with regions of apparent X-ray enhancement 
'clumps") and deficit (" bubbles"). The previously identifie d X-ray emitting "NE cloud" 



(~3" NE of the nucleus; lYoung et al.l 1200 ll : lOgle et al.l l2003l ) is clearly resolved. Here we 
focus on the inner ~200 pc (3") scale region that has been poorly studied in the X-rays 
before and now becomes accessible with the HRC deconvolved image (Figure Wp)- The 
subpixel resolution ACTS image of 0.1 s frame time was not used for deconvolution, because 
the appropriate model for the energy dependent, subpixel resolution PSF is not yet available 
and there is still low level (~5%) pileup in the nucleus (3.32 ± 0.048 counts s^^). 

The X-ray emission in the nuclear region is resolved into distinct components in the HRC 
image, namely a bright nucleus with elongation to the north, and two blobs ~0."7 from the 
nucleus at P.A.~6° and P.A.~54°, respectively. A bright knot of X-ray emission is uncovered 
at 1.4" to the NE of the nucleus at P.A.~25° (labelled as "X-ray knot N" in Figure [T]d, and 
throughout the paper referred as "HST-G" following HST features in the optical and avoid 
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inventing new nomenclature, see next section). A similar but fainter feature is located at 
1.3" SW of the nucleus (X-ray knot S, or "HST-H"). These two knots are also clearly present 
in the ACIS subpixel image (Figured]!), which shares a great similarity of morphology with 
the HRC images. 



3.2. Comparison with High Resolution [OIII] Image and Radio Maps 



We further compare the rich structures in [OIII]A5007 emission line and radio images 
to the X-ray emission. Figures and show the HRC image ove rlaid with coiitours o f the 
[OIII] emission from the ilf^T WFPC2/F502N image presented in ICapetti et all ( IiQQtI ) and 
the innermost zoom covers 150 pc regi on (2 arcsec a cross) . The labelled features following 
the HST clouds naming convention in lEvans et al.l (jl99ll ). Previ ous ACIS study al ready 
showed a general spatial correlation between the two wavebands ( lYoung et al.ll200ll ). the 
HRC image firmly demonstrates that there is remarkable morphological agreement between 
the X-ray clumps and the [OIII]A5007 emission, even at the smallest spatial scales. 



A similar compar ison with the VLA A-array 6 cm map (IWilson &: Ulvestad 



Gallimore et al.lll996al ) is shown in Figure[2t. We find that the X-ray "NE cloud" ( lOgle et al. 



1983 



20031 ) does not correspond to any enhancement in the radio lobe, although it is slightly offset 
(1" to the SW) from a bright radio filament. The newly identified X-ray knot (HST-G) is 
0.2" from the peak of a radio knot, where the coUimated radio outflow becomes more diffuse 
and lobe-like. 

The zoom (Figure |2]i) shows fine details of the spatial relations with the MERLIN 6 cm 
features (components NE, C, SI, and S2 in iGallimore et al.lll996al ). The X-ray elongation 
northward of the nucleus follows the linear radio emission between SI and C, then become 
faint and spatially anti-corr elated with the radio emission after the jet shows significant 
"bending" at component C ( IGallimore et al.lll996al ). 



3.3. X-ray Spectral Modeling of HST Cloud G 

The 0.1 s frame time ACIS image with subpixel resolution clearly resolves HST-G and 
HST-H (Figure [T]i), which are free of pileup. Because of the short exposure (1.4 ks), only the 
brighter HST-G have enough counts (386 ± 19) to attempt spectral fit using photoionization 
models. The X-ray spectrum was extracted using CIAO tool specextract from a 0.5"-radius 
circle centered on the X-ray knot, which was grouped to have a minimum of 1 5 counts per 



bin to allow usage of statistics. Spectral modeling was performed following IWang et al. 
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( l2011bl ) wit h the XSPEC p ackage (version 12.7.1) and photoionization model grids generated 
by Cloudy flFerlandl 12004 version 08.01). A single ionization (loglJ = -l.l^logNu = 21.7) 
component is not adequate to obtain a reasonable fit to the spectrum (reduced = 2-2). 
The addition of an optically thin thermal plasma component (APEC) with the best fit 
temperature of kT = l.OSto'il keV significantly improves the fit (xl = l-l; = 22; F- 
statistic=9.94 or probability of improvement 99.8%). There is a residual line feature at ~ 0.65 
keV (see Figure [3]) that is probably due to OVIII Lya, and can be fit with an additional high 
ionization component {logU = 0.9, xl = 0.8). The probability of improvement over current 
model is 90%. However, this has little effect on the thermal component. The flux contributed 
by the thermal component is -fo.5_2keV = 1.8 ±0.6 x 10~^^ erg s~^ cm~^ (-Z^o.5-2keV = 4.5 x 10^^ 
erg s~^), which is approximately 30% of the total X-ray emission from cloud G. 



Discussion 



The unprecedented high spatial resolution HRC image of NGC 1068 has enabled us to 
compare the sub-arcsecond structures of the X-ray emission and those of the HST NLR clouds 
in the nuclear region, where the ACIS data is affected by pileup. Overall, there is a good 
correlation between enhancements in [OIII] and X-ray emission (Figure [2^,b) in the nuclear 
region of NGC 1068, where the grating spectroscopic observations resolve the soft X-ray 
emission into many lines characterized by photoionization and photoexcitation (e.g., NVII, 
OVII, OVIII, and NelX), consis tent with previous suggestions that both originate from ga s 
photoionized by the AGN (e.g., IVoung et allboOll : ISrinkman et al.lbood : lOde et al.lbooah . 
This general morphological agreement appear s to be common for kpc-scale soft X-ray and 
[OIII] emission in nearby Seyfert galaxies (e.g.. iBianchi et al.ll2006U2010l : IWang et al.ll2011al : 
Paggi et al.ll2012l ). The ra tios between [QUI] and soft X-rays were found to effectively trace 



the ionization states (e.g., iBianchi et al.l 120061 : IWang et al.ll2009a[ ). Following these studies. 



we use [OIII] /soft X-ray ratios of different clouds to locate sites of shocks produced in strong 
jet-cloud interaction, as indicated b y low [OIII] /soft X-ray ratios due to enhanced thermal 
X-ray emission from shock heating (IBicknell et al.lll998l : IWang et al.ll2009al . |2011b|). 



Using the calibrated WFPC2 F502N image dCapetti et al.lll997h . we measured the [OIII] 
fluxes for the clouds (Evans et al. 1991) following WFPC2 narrow-band photometry in 
WFPC2 Data HandboolJ^. The 0.5-2 keV X-ray fluxes were measured from counts extracted 
from the same regions using the HRC image. The central positions (right ascension and 
declination in J2000) of these clouds and the rectangular regions used to estimate the [OIII] 



^,http : //documents . stsci . edu/hst/wf pc2/docmnents/heLndbooks/dhb/wf pc2_cover . html 
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and X-ray fluxes were listed in Table [T] It should be noted that there are finer [OIII] 
features resolved in the HST image that are not resolved in the X-rays, so these fiuxes only 
refiect the average values for the ionized gas. We assumed the phenomenological model 
(two-component bremsstrahlung plus power law) for the X-ray emission {kT = 0.45 keV and 
r = 1.0, Nh = 2.99 X lO^o cm-2; Young et al. 2001) and obtained Fo.5-2keV = 5.2 x 10^^^ 
erg s~^ cm~^ per HRC count s~^ for the fiux conversion. In Figure H] we show the [OIII] to 
soft X-ray ratio for five distinct cloud features along the jet (Tabled]) at various radii to the 
nucleus (~20 pc-100 pc). 

We note that three regions (HST-B+C, HST-G, and the previously unlabeled HST-H) 
show significantly lower [OIII]/X-ray ratio (~1) than the t ypical range spanne d by pho- 
toionized clouds (~3-ll in the Seyfert galaxies observed by iBianchi et al.ll2006l ). implying 



higher collisional ionization at these locations. We suggest that these clouds are strongly 
interacting with the jet, and that shocks driven into the obstructing clouds produce thermal 
X-ray emission. Indeed, our spectral fitting of cloud G revealed presence of kT = 1 keV 
hot plasma. The remaining two regions (HST-D+E and HST-F), which are clouds that 
bracket the MERLIN knot "NE" dCallimore et al.lll996b[ Figure Hi), show no evidence of 
jet-cloud collision. This scenario is consistent with the multi-wavelength morphology of the 
circum- nuclear region. A t the radio knot C, the P . A. of the jet changes from ~ 11° to ~ 33° 
( IGallimore et al.lll996a| ): iGallimore et al.l ( Il996bl ) suggested that it likely results from shock 
interaction between the radio i et and a dense molecular cl oud, which was indeed identified 
in 2.12/im H2(l-0) SI emission (jMiiller Sanchez et al.ll2009l ). Similarly, at the clouds HST-G 
and HST-H, the jet changes into a lobe after driving shocks into these clouds, as shown by 
our observations. There are other regions also show low [OIII] /X-ray ratios, such as the 
region between HST-D and HST-G, and a narrow "bar" feature (NW-SE) at the base of the 
radio lobe. These regions are likely associated with collisional gas as well, although much 
fainter and diffuse to measure the [OIII] emission. 

The kinematics of the ionized ga s also suggest jet - cloud impact. The long slit spec- 
troscopy of NGC 1068 with HST FOG dAxon et al.lll998[ ) found that the NLR gas at HST-G 
is strongly kinematica l ly per turbed and shows velocity systems separated by 1500 km s~^. 
Grenshaw fc Kraemerl ( 120001 ) show evidence for acceleration of the radio outflow from the 
nucleus and subsequently deceleration to systemic veloci ty. The same behavior was seen in 
the high ionizing potential coronal line region (GLR) gas ( jMiiller- Sanchez et al.ll201l[ ). where 
velocity tomography of coronal lines like [Si VI], suggested a turnover point for the radial 
acceleration at 80 pc NE from the nucleus, at the position of cloud HST-G, where the jet 
changes into a lo be. Moreover, the detec t ion o f extended highly ionized gas (e.g., [Fe VII]) 
in NGG 1068 led iRodriguez-Ardila et al.l (|2006[ ) to the suggestion that shocks related to the 
outflow could be required to power the line formation. 
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Last but not least, the spectral modeling of cloud HST-G revealed the presence of 
kT ~ 1 keV hot gas. Assuming a spherical symmetry for the gas (r = 0. 5^^ or 36 pc) , we ca n 



derive the emission measure and electron density for the hot gas (e.g., IWang et al.l l2009bl ) . 



Of most interest, the estimated thermal energy content deposited by the jet-ISM interaction 



1.9 X 1 ergs. Taking the velocity of the fast shocks driven by the jet to 



is Efj^ermal 

be y ~ 700 km s~^ (lAxon et al.lll998l ). the relevant timescale is the crossing time for the 



shocks to move through the 1" (72 pc) extent of the cloud, tc 



d/V = 10^ yr. If the 



jet-cloud interaction converts kinematic energy into heating of the hot gas, this implies a 
kinematic luminosity (or energy injection rate) of Lk.e. = 6 x 10^* erg s~^, which is lower 
but comparable to the kinematic luminosity measured from the jet-cloud interaction in NGC 
4151 {Lk.e. = 1 x 10^^ erg s~^). To put this in the context of the available jet power, we 
derive the total jet power of N GC 1068, -P,^ ^ ^. ^0^^ s~^, adopting the flux density 
of 4.8 ± 0.2 Jy at 1.465 GHz Jwilson fc UlvestadI EgSgf) arid the P... -P..^.. relation from 
a sample of radio loud galaxies with X-ray cavities (jCavagnolo et al.ll2010l ). A negligible 
amount of the jet power {Lk.e./ Pjet = 10~^) is lost to the ISM during the encounter with 
the cloud G. 

Our Chandra HRC and ACIS results show three locations of intense interaction between 
the radio jet in NGC 1068 and the NLR clouds in the inner 20 pc: the jet first encou nters a 
dense clump at knot C (HST-B+C), produces a shock front ( Gallimore et al. 1996a| jbl) and 
changes direction; then it sweeps through between clouds D-fE and F with no evidence of 
interaction; it interacts again with clo uds HST-G in the N E and HST-H in the SW, changing 
morphology and ionization structure (ICapetti et al.lll997l ). 



5. Conclusions 

We have obtained high spatial resolution X-ray image of the nucleus of NGC 1068 using 
the HRC, which allows, for the first time, a direct view of the innermost 1 arcsec (72 pc) 
radius region. 

The HRC image resolves the narrow line region into X-ray emission clumps matching 
bright emission-line clouds in the HST [OIII] A5007 images. Two distinct X-ray knots are 
revealed at 1.4 arcsec NE and SW of the nucleus. With image deconvolution technique, we 
resolve sub-arcsec X-ray emission features tracing the coUimation and bending of the jet. 

Based on the combined X-ray, [OIII], and radio continuum morphology, we identify 
three locations of radio jet-cloud interaction. The [OIII] to soft X-ray ratios show that these 
clouds are strongly affected by shock heating, whereas near two clouds the jet simply thrust 
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through with no signs of strong interaction. This is further strengthened by the presence of a 
kT ~ 1 keV coUisionally ionized component in the ACIS spectrum of the cloud G. We further 
estimate that the kinematic luminosity (or energy injection rate) of the jet-driven shocks is 
Lk.e. = 6x 10'^^ erg s~^, which shows a negligible amount of the jet power {L^.e./ Pjet = 10"'^) 
is lost to the ISM during the jet-cloud interaction. 



Our results show that the NGC 1068 radio jet plays a crucial role in shaping the morphol- 
ogy, the kinematics and the ionization structure of the innermost NLR. Recent high spatial 
resolution studies based on integral-field spectroscopy with Adaptive Optics in the near-IR 
start to provide velocity fiel ds and trace highly ionized gas in the radial outflows in the 



NLR of S eyfert galaxies (e.g..lStorchi-Bergmann et al.ll2010l : iMiiller-Sanchez et al.ll201l[ ). As 



shown in IStorchi-Bergmann et al.l ( 120101 ). [Fell] may be enhanced where jet-induced shocks 
probably release the Fe locked in grains. We expect to see enhanced [Fell] peaks at the 
locations of jet-NLR interaction suggested in this work. In addition, the deep ACIS HETG 
grating observations of NGC 1068 may provide a more detailed picture of the kinematics 
and ionization of the hot phase ISM in NGC 1068. 
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Table 1. Measured X-ray and [OIII] Fluxes 



Cloud 
label 


Distance 
toNuc. (") 


Distance 
to Nuc. (pc) 


R.A.=02h 42™ 
(J2000) 


Decl.=-00° DO' 
(J2000) 


Extraction 

Area 


[OIII] flux 
(c.g.s) 


0.5-2 keV flux 
(c.g.s) 


[OIII] /soft X 


B+C 


0.25 


18 


40.=719 


47''71 


0" 4x0" 45 


15.4 


16.6 


0.9 dz 0.1 


D+E 


0.6 


43 


40.=724 


47^' 18 


0'.'4x0'.'4 


14.2 


1.9 


7.3 ±0.6 


F 


0.7 


50 


40.=752 


47^' 40 


0'.'35x0^'35 


14.0 


1.5 


9.6 ± 1.0 


G 


1.4 


101 


40.^758 


46'.'59 


0'.'45x0'.'4 


7.7 


2.9 


2.6 ±0.2 


H 


1.3 


94 


40.^675 


49'.'31 


0'.'45x0'.'4 


0.73 


0.38 


1.9 ±0.9 



Note. — The [OIII] cloud designations are adopted from Evans et al. (1991), except cloud H was not previously labelled and assigned here. A uncertainty of ~5% for [OIII] 
flux is adopted (Capetti et al. 1997). The [OIII] and X-ray fluxes are in unit of 10~^^ erg cm~^ per arcsec^. 
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ACIS-S native pixel 
1" 
72 pc 





Fig. 1. — (a) The HRC image of the 7 arcsec across (500 pc) nuclear region of NGC 1068; (b) EMC2 
deconvolved HRC image; (c) The ACIS 0.3-2 keV image (0.1 s frame time) of the same region, displayed at 
the native pixel (0.492"); (d) The subpixel repositioned, rebinned ACIS image of the same region displayed 
at 1/4 of the native pixel (0.123"). A 1" scale bar is shown in all panels, corresponding to 72 pc at the 
distance of NGC 1068. 
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Fig. 2. — (a) The deconvolved HRC image of the X-ray emission in NGC 1068 overlaid with the contours 
of HST/WFPC2 [OIII] A5007 hne emission (Capetti et al. 1997). Contour levels are from 3 x 10'^^ to 
3 X 10^^^ crs cm^^ s^^ (PC pixel)^^. (b) Zoom-in image showing the detailed correspondance of [OIII] and 
X-ray features in the innermost region; (c) The deconvolved HRC image overlaid with the contours of VLA 

5 GHz emission (Wilson k Ulvestad 1983). The beam size is 0'.'49x0'.'38. Contour levels are 0.3, 0.7, 1.4, 
2.8, 5.6, 16, 27.9, 59.5, 120, and 224 mJy beam'^ (d) The innermost region overlaid with the MERLIN 

6 cm contours (Gallimore et al. 1996). The beam size is 0'.'065. Contour levels are 0.862, 1.73, 3.50, 7.06, 
14.23, 28.7, and 57.8 mJy beam~^. Note the 0'.'5 scale bar in the zoom-in panels (b) and (d). 
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Fig. 3. — The X-ray spectrum of the cloud HST-G, and the best fit model consisting of a 
kT =1.08 keV thermal component (the dash-dotted line) and a low ionization photoionized 
emission component (the dotted line). 
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Fig. 4. — The [OIII] to soft X-ray ratio as a function of the cloud's distance to the nucleus. 
The (red) diamonds are the clouds as labeled in Figure [2] and Table [1], and previous mea- 
surements for NGC 4151 clouds (Wang et al. 2009a; Wang et al. 2011b) are also shown 
for reference. The shaded area is the range of ratios in Bianchi et al. (2006) that can be 
explained with photoionization modeling. The blue, gold, green, and magenta lines indi- 
cate the [OIII]/X-ray ratios for the NLRs of NGC 7212, NGC 3393, Mrk 3, and NGC 1386 
(Bianchi et al. 2006), respectively. 



